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ABSTRACT 


An experiment was performed to empirically determine the fading characteristics 
of backscattered radar signals from four agricultural targets at 9 GHz. After a short 
review of the statistics of Rayleigh fading backscatter, the data processing method and 
results of the data are analyzed. Comparison with theory shows adequate agreement with 
the experimental results, provided of course, the targets are modeled in a correct 
manner . 
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1.0 INTRODUCTION 


In studying the scattering of radar signals, fading is seen to play an important 
role. Fading, in the case of any coherent signal , is caused by the constructive and 
destructive interference of the signal components reflected by individual scatterers 
within a resolution cell . Thus, the phenomena of fading will be a function directly 
dependent on the target being studied. 

This fading, or fluctuation, observed in the radar return introduces a significant 
complication in studying the scattering properties of a particular target. Making a 
single instantaneous measurement of scattered power will introduce errors since there 
are no means by which we can determine how much fading is indeed taking place. Thus 
it is necessary to make a number of measurements which can then be averaged so that 
the effects of fading are reduced. 

The phenomena of fading appears in most radar studies and its effect produces 
a certain amount of uncertainty in the measurements being made. In the case of an 
imaging radar, the effect of fading is to cause speckle in the image. But it has been 
shown that even when very little averaging is performed by the imaging system, the 
eye-brain system performs an adequate amount of averaging when viewing a radar 
image [1]. In detection radar systems, the problem of determining the probability of 
detecting a fading target has been studied extensively by Swerling[2j and a variety of 
other investigators [3,4}. 

Because scatterometers are by definition calibrated devices used in determining 
scattering information on an absolute scale, fading poses particular problems. The 
question that must be asked is: How many samples of the scattered power must be averaged 
before certain confidence limits can be placed on the final average? This is often a 
difficult question to answer unless some a priori knowledge of the scattering properties 
of the target is known. This knowledge may be either theoretically derived or it may be 
based on data previously collected, de Loor [5] made an attempt to characterize the 
fading properties of certain agricultural targets although his experiments are limited 
to incidence angles larger than 80°. 

In order to shed some experimental light on the subject of fading, it was 
decided to conduct an experiment to determine the fading properties of four selected 
agricultural targets. This report presents the results of this experiment. 



2.0 EXPERIMENTAL INSTRUMENTATION AND TECHNIQUE 

The radar used in Hi is study was a wide band 8-18 GHz radar spectrometer 
scatterometer [6]. This is an FM-CW system employing two antennas (transmit and 
receive) mounted atop a 26 meter hydraulically operated boom. The boom is in 
turn mounted on a truck for mobility. The transmitted signal is frequency modulated 
by means of a triangular waveform with an amplitude adjusted to provide a frequency 
deviation of 400 MHz. Table 1 presents a summary of the basic system specifications 
with a system block diagram shown in Figure 1 . 

TABLE 1. 

BASIC RADAR SPECIFICATIONS 


TYPE 

FM-CW 

Modulating Waveform 

Triangular 

Frequency 

8-18 GHz 

FM sweep: Af 

400 MHz 

Transmitter Power 

10 dBm (10 mW) 

Intermediate Frequency 

60 kHz 

IF Bandwidth 

3.58 kHz 

Antennas 


Height above ground 

26 m 

Reflector Diameter 

61 cm 

Feeds 

Cavity backed. 


log perodic 


Effective Beamwidths of Product Patterns 
(degrees) 

Azimuth Elevation 

2.94 3.43 

The system was calibrated against a Luneberg lens of known cross section. 

Because of technical difficulties, fading measurements were made only at 9 GHz with 
W polarization. 


equency Calculated Antenna Gain 

■(GHz) (as) 

9 31 .2 
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FIGURE 1. BASIC BLOCK DIAGRAM OF THE 8-18 RADAR SPECTROMETER. 






The measurements were made by locating the trucks at the edge of the field 
containing the target of interest. About 140 measurements were collected from each 
field at each of six incidence angles: 0°, 10°, 20°, 30°, 50° and 70°. After each 
individual measurement the boom was moved a few degrees in azimuth in order to 
view an independent resolution cell . Non-overlapping cells were assumed independent. 
Targets studied were corn, alfalfa, soybeans and bare soil. 


3.0 AVERAGING PROPERTIES OF PANCHROMATIC RADAR 


The term "panchromatic", although usually applied to a multi-frequency signal 
in the visible portion of the spectrum can just as well be applied to a broad spectrum 
radar signal. Since the 8—18 GHz radar used in the experiments under discussion is 
panchromatic it may be helpful to review some of the averaging properties of a 
panchromatic radar. 

In 1963 Birkemeir and Wallace [7] considered the case of a target whose 
elements are uniformly distributed along a line of length L as shown in Figure 2. 

After deriving an expression for the frequency autocorrelation function, they chose 
the first zero as the point of decorrelation. The corresponding frequency spacing for 
independence was found to be given by: 

i,= rH? MH * o°> 

where 0 is the angle of incidence. 

Other variations of this problem have been studied by Ray [8], Moore, et al. 
[9] and Waite [10] and in all cases the results yield expressions similar to Eq, la. 

Ray's [8] model is somewhat more applicable to vegetation because his target is a 
volume having a maximum dimension D (as measured radially from the antenna); his 
result takes the form: 

Af= MHz (lb) 
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4.0 A THEORETICAL DESCRIPTION OF RADAR BACKSCATTER 


Before proceeding to a discussion of the experimental results, it will be 
helpful if we review some of the basic statistical properties of radar backscatter. 
Following a development provided by Beckmann [1 1] let us consider the returned 
E field, E s , from a single scatterer: 






- 




V 


( 2 ) 

(3) 


where 




(4) 


Now let us assume E and / 
between ip and \J> + 2 tt . 

independent scatterers: 


to be random variables and that </ is uniformly distributed 
Consider the total E field returned from a collection of 


•=T- zr y s*£ 
* 2 - ft <2 


(5) 
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where n is the number of scatterers illuminated, R is the magnitude of with angle 
0. The illustration in Figure 3 may be helpful. 



Figure 3. 


We can easily see that the following is true: 


1 j*. 


Y= £ & •=• 



(6a) 

(6b) 


Now, if we assume independence of X^'s and *s and if we let n become large we 

can apply the central limit theorem [12] and thus conclude that x and y are normally 
distributed with means: 

<*> = > 

(7a) 

<Y ' > ~ \z-, <4. > 

(7b) 

Also, if the E^ and </ are uncorrelated we can write 



(8a) 
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Z7T 

-£r/ (8b) 
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Similarly, 

< Y> 

Likewise we can calculate the variances 


x - <x'> *ZT <eY><co^d> = (9a) 

/ * / 


■=7' i = </*> *21 i> ^z’ - <^7 (5b) 

£ -/ 

Thus, in accordance with the central limit theorem we can write expressions for the 
probability density functions of X and Y as follows: 


/ 

— X^/z crj- * 


' y Z7T' 


(10a) 

/ 




<^ / 7 

(10b) 


Now let’s check to see if X and Y are correlated : 

*>= 2f: 2ZT acs d <*, w 6 > <") 

i. ~ / -r ■= * ' * 

By our assertion that the E's and /'s are uncorrelated we can write 

0 2 ) 

since cos sin >=0 . 
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Thus, 
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fs Vy z ) / /z c/~ *- 


^77" C7 — 

Now, making the transformation to polar coordinates [13] we see that 


— - — 5 — 
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where 
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Note that 
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Thus, 




-**£ 


and 


a77- 


Jo 


211/ ~ 


which is a Rayleigh probability density function. 

Next we desire to determine the distribution of the return power, 
only a matter of making an appropriate change of variables : 

w-= 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

This is 

( 20 ) 
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where W is the return power, V is the envelope voltage amplitude of the received 
signal and K is a proportionality constant. Dropping the constant K for convenience 
and making the change of variables it can again be shown [13] that: 

<2,) 

Although Eq. 21 is in reality an exponential density function it is often referred to 
as a Rayleigh power density function for obvious reasons. 

On closer observation one will see that Eq. 21 is also a chi-square 
density function with two degrees of freedom. The chi-square density function is 
defined as : 



where n is the number of degrees of freedom. The mean and the variance of the chi" 
square distribution are given as; 

— C7— ' * = ^ (23) 


where a is the standard deviation of the Rayleigh power density function. 

Next let us consider the average of N samples, each distributed, as we 
have shown, with a chi-square distribution with two degrees of freedom. We then 
have 




(24) 


Consider 


23 (*'/' 
/ = / 


(25) 


By the addition theorem of the chi-square distribution [14] we know the distribution of 
Eq. 25 will be chi-square distributed with 2N degrees of freedom. Thus, if n = 2N 
we can write the density of Eq. 25 as: 



/ 








r 


(26) 


9 



What we actually wish to determine, however, is the density function of W where 


IV - 



(27) 


Making a change of variables we see 


7 #/ — - /I C 


^ C K 
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A/ VJ I , — r A/ - v 




(28) 


which is a Gamma density function for W. 

Finding the expected value and variance of W in the usual fashion we obtain: 


^Trz7- 

/£_ <2 cr~ 
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(29) 

vfaj* 

W /l/ 

(30) 


5.0 DETERMINATION OF THE NUMBER OF INDEPENDENT SAMPLES AVERAGED 


At this point we have shown that we expect the distribution of the return power 
from a target consisting of a complex of discrete scatterers to be describable by a 

Gamma density function. We will follow this assumption throughout the remainder of this 
report. 


We must now determine a method from which the actual number of independent 
samples, N,can be calculated. From Eqs. 29 and 30 we see how this can be done. That is. 


C&jdl _ 




— 4/ 


Vv 


(31) 


Thus, we can determine N, the number of independent samples averaged if we can 
determine and cr-^^ . These parameters were estimated in the usual fashion 

by employing Eqs. 32 and 33: 


Wr-j - 
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21 
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(32) 
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where m is the number of measurements taken. 

In the cases of corn, soybeans, and alfalfa, the number of measurements averaged 
was 140 and for bare ground 152 samples were averaged. 

If we follow Eqs, 31 through 33 we can determine N. The results are 
shown in Table 2. 


TABLE 2. 

CALCULATED NUMBER OF INDEPENDENT SAMPLES, N, 
BASED ON MEASURED DATA 


Incidence 

Angle 

0° 

10° 

o 

o 

CN 

o 

o 

CO 

50° 

70° 

Bare Soil 

2.5 

1.3 

2.3 

3.7 

3.8 

7.8 

Com 

1.0 

1.3 

7.0 

6.0 

8.5 

11.3 

Soybeans 

2.8 

2.7 

2.2 

5.0 

4.2 

10.0 

Alfalfa 

1.7 

2.1 

1.8 

2.5 

5.2 

8.3 


At the time of measurement, soil samples were taken in three layers to 
determine moisture content by weight. These layers as measured from the soil 
surface are: 0-1 cm, 1-2 cm and 2-5 cm. The per cent moisture is shown in Table 
3. 


TABLE 3. 

SOIL MOISTURE CONTENT BY WEIGHT 


Depth (cm) 

0-1 

1-2 

2-5 

Bare Soil 

19% 

20% 

21% 

Corn 

26% 

24% 

25% 

Soybeans 

6% 

11% 

26% 

Alfalfa 

18% 

25% 

19% 


11 


Let us now calculate the number of independent samples as predicted by Eq. 

1b. The range resolution of the scatterometer used in this experiment is given 
in meters by • 
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(P 


(37) 


where © is the incidence angle. Thus if the maximum dimension, D', of the target 
as measured radially from the antennas is less than AR, we will assume D = D 1 
where D is defined in Section 3. If, however, D' is greater than AR we will 
assume D = AR. If complete penetration through the vegetation is assumed, then 
D' = Hsec0, where H is the vegetation height. For cases where no penetration is 
expected, D' is the difference between the maximum and minimum radial distances to the 
target as determined by the antenna beamwidth. Shown in Table 4 are the heights 
of the various crops at the time of the experiment. 


TABLE 4. 


Crop 

Crop Height 

Soybeans 

0.36 m 

Corn 

2.45 m 

Alfalfa 

0.63 m 


With the information of Table 4 and employing Eq. 37 it is possible to determine 
Af as defined by Eq. lb. Ray's model [8] was chosen as it seems to be 
representative of the models discussed in Section 3.0. Knowing Af we can 
determine N by dividing the modulation bandwidth, 400 MHz, by the frequency 
spacing between independent samples, Af. The results of this calculation are shown 
in Table 5. It was necessary to make certain assumptions concerning the amount 
of penetration taking place. Since corn and soybeans are both row crops, the radar 
return will normally have contributions from scatterers at various depths between the 
top of the plants and the ground. Hence total penetration was assumed. Alfalfa, 
however, is not a row crop and it was not immediately obvious as to what degree 
penetration occurred. For this case Table 5 shows two entries for alfalfa. The first 
entry, N z , assumes zero penetration while the second, N^, assumes total penetration. 
In the case of bare ground, no penetration was assumed due to its small skin depth 
at 20 per cent moisture content [15] . 
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TABLE 5. 

NUMBER OF INDEPENDENT SAMPLES, 
N, AS CALCULATED USING EQ.lb 


Look 

Angle 

Corn 

N 

Soybeans 

N 

Bare Ground 
N 

Alfalfa 

N 

z 

N , 

0° 

4.08 

0.959 

0.026 

0.026 

1.7 

10° 

4.15 

1.6 

0.694 

0.694 

2.3 

O 

o 

<N 

4.34 

2.42 

1.39 

1.39 

3.2 

30° 

4.72 

3.50 

2.08 

2.08 

4.4 

40° 

5.34 

5.25 

4.16 

4.16 

6.1 

50° 

6.34 

6.34 

6.34 

6.34 

6.34 

o 

o 

o 

8.16 

8.16 

8.16 

8.16 

8.16 

0 

o 

11.76 

11.76 

11.76 

11.76 

11.76 


The reason for N being the same for all crops at angles past 40° is that the range 
resolution of the radar is dependent solely on incidence angle for angles of 50°, 60 
and 70°. 


6.0 DATA INTERPRETATION 


In order to facilitate interpretation of the data, it is presented in a graphical 
format of three types. The first is in the form of histograms of the measured scattering 
coefficient, ct°, and are presented in Appendix A. To more easily see the 
effects of averaging, each point was divided by the mean value of o at each angle 
for that particular target. This insured that the mean occurred at a value a° = 1 
for all histograms. Although this is merely a qualitative tool, it does show the effect 
of averaging as the incidence angle increases from 0° to 70°. Note that the 
distributions do indeed look chi-square in appearance, as was predicted in Section 
4.0. Figure 4 shows the mean value of a° plotted versus incidence angle for each 
of the four targets under discussion. 
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o° in dB p , o° in dB 



Crop Type: Corn 
Frequency (GHz): 9.0 
Polarization: W 
% Soil Moisture by Weight: 26.0 

Crop Height (cm): 245 



W”20 30 40 50 60 70 

Incidence Angle in Degrees 



Crop Type: Soybeans 
Frequency (GHz): 9.0 
Polarization: VV 
% Soil Moisture by Weight: 6 
Crop Height (cm): 36 
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Crop Type: Alfalfa 
Frequency (GHz): 9.0 
Polarization: VV 
% Soil Moisture by Weight: 18. 0 
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Crop Type: Bare Soil 
Frequency (GHz): 9.0 
Polarization: VV 


% Soil Moisture by Weight: 20.0 
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Figure 4. 


Scattering coefficient, a°, versus incidence angle, 9. 


a) corn b) alfalfa c) soybeans d) bare soil 



The last melTiod of dafa presentation which gives a more quantitative 
description of the averaging process is that of plotting the data in Tables 2 and 5 
on the same scale. This is done in Figures 5 through 8, Since the histograms do 
have some empty spaces in them, it was felt that a smooth curve would be a better 
representation than actually joining the points by lines in a rather absolute manner. 

It should be noted that when Eq. lb predicts less than one sample, the value was 
elevated to a value of one since less than one sample is physically unreal. 

We begin by noting that Ray's model [8] is a fairly good predictor of N for 
soybeans and bare ground. The trend is very close to the trend of the experimental 
results and the absolute values of the two curves are fairly consistent with one 
another. This is not the case for corn, however. The experimental curve for corn 
exhibits a sharp dip for values of 0 less than 25°. This is attributed to an incorrect 
assumption concerning the degree to which the corn contributes to the measured return. 
If, for example, we assume that at small incidence angles the return is due mainly to 
the soil, we can understand why the curve exhibits the tendencies we see. This is to 
say that at small incidence angles we are "seeing" mostly the underlying soil. At 
larger incidence angles, however, the corn significantly attenuates the signal so that 
the effect of the soil is insignificant to the total return. 

This hypothesis is supported by two facts. First we see in Table 3 that the soil 
was very wet underneath the com; 26 per cent moisture by weight. This fact, combined 
with the fact that corn is a row crop and thus not as dense as a crop such as alfalfa or 
wheat, tends to augment the effect the soil will have on the scattered signal. Secondly, 
we can see from Figure 9 the effect of increasing the incidence angle on the return from 
corn with wet underlying soil . These curves show the tendencies of cr° as a function of 
look angle for soils with moisture contents of 42 .0%, 18.7% and 4.5% [17], The 
tendency for the curves to merge at angles greater than 30 indicates the signal from the 
soil is greatly reduced by the corn at these higher angles while at smaller angles the soil 
return dominates. 

Looking at the experimental and theoretical curves for alfalfa we see another 
interesting result. Note that two curves for Ray's model have been plotted. One 
assumes no penetration into the crop while the second assumes penetration down to 
the ground. At zero degrees, both curves predict the same result and they also agree 
with the experimental results. As the incidence angle increases, however, the curve 
assuming total penetration begins to deviate quite drastically from the experimental 
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curve while the second theoretical curve tends to be comparable to the measured 
result . 

Thus it seems that scatter from alfalfa tends to be a result of surface scatter 
rather than volume scatter. This effect had been noted earlier by Ulaby [16] in 
his studies of the angular dependence of ar° for different crop types. 


7.0 CONCLUDING REMARKS 


An experiment was performed to determine the fading characteristics of radar 
backscatter from agricultural targets. The assumption was made throughout the study 
that the distribution of the fading amplitude of radar backscatter from vegetation can 
be closely approximated by a Rayleigh distribution. For this assumption to be correct 
the target must consist of many random and independent scatterers. For naturally 
occurring vegetation such as a wooded area this is probably a good assumption. For 
agricultural vegetation, however, this may not be the case. Corn, soybeans and even 
bare cultivated fields are normally prepared and planted in a very orderly fashion by 
the farmer. The resulting geometric regularity produced may then make the validity of the 
Rayleigh assumption somewhat questionable. Schuchman and Drake [18] have reported 
that row orientation plays an Important role in radar return from corn. The assumption 
may also break down for bare cultivated soil in which case a good probability of 
specular scatter exists at angles near nadir. 

Another phenomena that should be investigated is that of signal refraction by 
the vegetation. This may be particularly important at angles away from nadir in which 
case a small amount of refraction may produce a relatively large change in the effective 
target extent. In addition to refraction, target penetration has been shown to play a 
significant role in altering effective target extent, thus affecting the fading process. 

Implementation of a non-coherent radar would provide more quantitative 
information on the range of validity of the Rayleigh assumption. A comparison of the 
fading process as measured by such a system to that as measured by a coherent system 
would then answer many of the questions as to the types of targets which have previously 
been assumed to scatter radar signals incoherently. 
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APPENDIX A 


HISTOGRAMS OF THE POWER DISTRIBUTIONS OF 
RADAR BACKSCATTER FROM BARE GROUND, SOYBEANS, CORN AND ALFALFA 
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Figures Ala.- Alb. Histograms of a 0 for corn normalized byo® the mean value of a° 
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Figures Ale. - Aid. Histograms ofo°for corn normalized bya° the mean value ofo° 
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Figures Ale. “Alf. Histograms of<r°for corn normalized by a°, the mean value ofcr° 
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Figures A2a*~ A2b* Histograms of <j° for soybeans normal! 
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Figures A2c.“ A2d. Histograms of a° for soybeans normal! 
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